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The field-driven reorientation transition of an anisotropic ferromagnetic monolayer is studied 
within the context of a finite-temperature Green's function theory. The equilibrium state and the 
field dependence of the magnon energy gap Eq are calculated for static magnetic field H applied in 
plane along an easy or a hard axis. In the latter case, the in-plane reorientation of the magnetization 
is shown to be continuous at T = 0, in agreement with free spin wave theory, and discontinuous at 
finite temperature T > 0, in contrast with the prediction of mean field theory. The discontinuity 
in the orientation angle creates a jump in the magnon energy gap, and it is the reason why, for 
T > 0, the energy does not go to zero at the reorientation field. Above the Curie temperature Tc, 
the magnon energy gap Eo(H) vanishes for H = both in the easy and in the hard case. As H is 
increased, the gap is found to increase almost linearly with H , but with different slopes depending 
on the field orientation. In particular, the slope is smaller when H is along the hard axis. Such a 
magnetic anisotropy of the spin wave energies is shown to persist well above Tc (T as 1.2 Tc). 

PACS numbers: 75.70.Ak, 75.30.Gw, 76.50.+g 

I. INTRODUCTION 

A great experimental achievement is the ability to grow epitaxially well defined monolayer and sub-monolayer 
ferromagnetic films. Such films serve as model systems in which to study basic aspects of magnetic ordering in low 
dimensions as well as to test theories of magnetism. Results gleaned from these studies have wide impact across 
several areas from critical phenomena to electronic band structure theory. There are also important implications 
for applications because the issues involved lie also at the heart of problems in interface magnetism. A problem of 
particular importance is the formation of magnetic anisotropies in thin transition metal magnetic films. 

Critical phenomena in low dimensional magnets and the formation of magnetic anisotropies are linked together at 
a fundamental level. Long range order in two dimensional magnets depends upon the long wavelength behaviour of 
the Goldstone mode of the spin system. This behaviour is strongly affected by the existence of a zero momentum 
energy gap. In most real systems, the largest contribution to this gap comes from second order corrections to the 
exchange energy associated with electronic correlations creating the ordered state. These corrections appear in the 
form of magnetic anisotropies determined by spatial symmetries of the local atomic environment. 

Anisotropy formation is a higher order effect in the sense that the energies involved are much less than the energy 
associated with the ordering temperature Tc- In two dimensional systems, the local fields responsible for magnetic 
anisotropies are very different from the local fields existing in three dimensional bulk systems, and the same reduction 
in dimensionality has a profound effect on the critical behaviour of the system. In real magnetic systems accessible 
to experimental study, the formation of magnetic anisotropies and the critical behaviour of the magnetization are 
therefore inextricably linked. 

In this paper we show that there are striking features associated with magnetic anisotropies which appear in 
spectra produced by spin wave excitations in two dimensional ferromagnets. There are very clear and measurable 
finite temperature effects produced by anisotropies that can be observed below the ordering temperature, and we 
demonstrate that one can also observe easy and hard directionality in magnetic excitation energies above the ordering 
temperature. 

A few experiments in epitaxial ultrathin ferromagnetic films have explored the nature of the phase transition and 
the persistence of magnetic anisotropies well above Tcr^ a similar effect was observed in quasi- two-dimensional 
easy-plane ferromagnets^ The role of the energy gap in stabilizing ferromagnetic order in two dimensions is well 
known in critical phenomena, but less well understood in terms of experiment. Experimentally, the gap can be 
observed from the spectra of long wavelength spin waves, and can be most clearly seen by taking the magnetic system 
through a reorientational transition with the application of a static magnetic field along a hard anisotropy direction^ 
Theoretically, the equilibrium configuration can be determined, e.g. by the vanishing of the total static torque acting 
on spins: for small fields the magnetization orients between the direction of the field and the closest symmetry axis of 
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the anisotropy until, at some specific critical magnetic field strength H c , a reorientation of the magnetization along the 
field direction occurs. As regards the spin wave excitations away of the equilibrium state, at low temperatures their 
energy can be easily calculated, using the Holstein-Primakoff transformation from spin operators to boson operators^ 
in the framework of a non interacting spin wave theory^*2*SiS For H = the spin wave energy is essentially determined 
by the anisotropy but, as H increases, the energy barrier to deviate away from the easy axis gradually decreases until, 
for H greater than H c , the field itself constitutes an energy barrier. This results in a minimum value for the energy 
of a spin wave at H c . Within free spin wave theory, valid at low temperatures, this minimum is predicted to be 
exactly zero for zero wavevector^ It is the aim of the present paper to investigate, at higher temperatures, the 
field dependence of the equilibrium configuration and of the energy of the magnetic excitations of a ferromagnetic 
monolayer with a uniaxial anisotropy. Two different theoretical approaches have been used, namely mean field theory 
(MFT) and Green's function theory (GFT). The different predictions of these theories are exposed in Sections II and 
III, respectively, both above and below the Curie temperature Tc- Within MFT a constant magnetization modulus 
is predicted for field smaller than the reorientation field, whereas within GFT a more accurate treatment of spin 
fluctuations leads to a field dependence of the magnetization modulus in the same field range. As a consequence, for 
finite temperature < T < Tc, we find that the two theories give qualitatively different results as regards the order of 
the field-induced transition: within MFT the reorientation of the magnetization is continuous and the magnon energy 
gap vanishes at the reorientation field; within GFT the reorientation is discontinuous and the gap does not vanish at 
H c . Above Tc, both theories are able to account for anisotropy effects on the energy of the magnetic excitations, but 
quantitative differences are present. Finally, the conclusions are drawn in Section IV. 



II. MEAN FIELD THEORY 



We consider a magnetic monolayer with Hamiltonian 
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where J > is the nearest neighbor ferromagnetic exchange interaction, K 2 > is a uniaxial single-ion anisotropy 
favoring the Z direction in the monolayer plane (XZ), and H = (H x ,0, H z ) is an external magnetic field applied 
in plane: see Fig. 1. The last term in Eq. (JTJ is the magnetic dipole-dipole interaction, favoring the alignment of 
the magnetization within the film plane; denotes a vector joining two different lattice sites and a is the lattice 
constant. In the following we will assume a simple quadratic lattice, so that the number of nearest neighbors is 4. 
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A. Equilibrium condition 

The effective field H Bit acting on a spin in the ferromagnetic film is composed of an exchange field H ex 
a shape demagnetizing field of modulus 4Trg^B(S z ) favouring the film plane (Y = Q); a uniaxial anisotropy field 
Hk = 2A g 2 ^f ^ favouring the crystallographic Z direction in the film plane; an external field H applied within the 
film plane, along the easy (Z) or the hard (X) axis. By (S z ) and (S z ) we denote, respectively, the thermal average 
of the spin component referred to the crystallographic reference frame X, Y, Z and to the local reference frame x, y, z 
(see Fig. 1). More precisely, z is the equilibrium direction of the in-plane magnetization, forming an angle <f> with 
the easy in-plane crystallographic direction Z . In the mean field approximation, the thermal average (S z ) referred to 
the local frame can be calculated self consistently in terms of the classical Langevin function (S z ) = SC(x), where 
C(x) = cothz - i and x = g/iB , H °»" s ' . 

The exchange field is parallel to the magnetization while the shape demagnetizing field does not provide any 
preference for a particular in-plane direction. The equilibrium value for the angle <j> that the magnetization forms 
with the crystallographic easy axis Z is determined by requiring the vanishing of the static torque M x (H + i/fj-Z). 

For field along the easy axis, H = H Z Z, there is no competition between the applied field and the anisotropy field, 
so that = for any field value. The thermal average of the magnetization is obtained solving (S z ) = SC(x) with 
x = x eas y = ^(g^H 2 + 2K 2 {S Z ) + U(S Z )). 
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the condition for vanishing static torque is 



For field along the hard axis, H = H X 'K, the modulus of the anisotropy field is Hk — — 1 ^ n ^ s case : 
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FIG. 1: Crystallographic (X , Y, Z) and local (x,y,z) reference frames used in this paper. The uniaxial in-plane anisotropy K2 
favors the Z direction within the film plane (Y = y = 0) while a magnetic field is applied in plane along the easy (Z) or the 
hard (X) direction; <f> denotes the angle that the magnetization M, directed along the local axis z, forms with the easy in-plane 
direction Z. 



The thermal average is obtained solving the self consistent equation (S z ) = SC(x) with x = x hard = 
-^ T (g^ B H x sm4> + 2K 2 {S z ) cos 2 <f> + 4J(S z )). As H x is increased from zero to the critical field value H x = 2K ^) ^ 
it turns out that <f> ranges continuously between and tt/2 and the magnetization modulus M — (S z ) is independent 
of the field intensity in the whole field range < H x < H x (see Fig. 2a). 

The latter feature can be readily accounted for observing that, upon substituting the equilibrium condition 

appropriate for < H x < H x (namely <j> — sin -1 ( I^L ) ) in the expression for x hard , one obtains that 
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the magnetization modulus M — (S z ) — S£(x) is determined by the solution of the self consistent equation 
x = Y^?p{2K2 + 4J)(cotha; — §), (where the field intensity H x does not explicitly appear) in the field range from 
H x = up to H = H x . In contrast, for H x > H x the magnetization is field dependent since one has <fi = tt/2 
and x — Y^[gfiBH x + 4JS'(cotha; — i)]. Notice that the field dependence of M for H x > H x is not apparent from 
Fig. 2a since it is very weak for the chosen parameter values. 

For zero field, one has x = (4 J + 2K 2 )(S Z ) and the Curie temperature Tc can be calculated, as usual, by 
requiring that the initial slope of C{x) is equal to 1, i.e. to the slope of (S z ) / S. One obtains 
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B. Resonance frequency 



The resonance frequency can now be calculated by solving the linearized Bloch torque equations dM/dt — — 7M x 
H aff where 7 is the gyromagnetic ratio. For field along the easy axis, H = H Z, the time varying magnetization can 
be assumed to be of the form M = ^Xmjf + Ymy^j e~ luJt + 7ig/i B (S z ), while the effective field is 
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We note that the effective exchange field was not included in Eq. J2J since it is parallel to the magnetization, thus it 
does not contribute to the uniform mode frequency. The resonance frequency is then readily found to be 
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FIG. 2: Comparison between the results of mean field theory (MFT, left) and Green's function theory (GFT, right) for the 
magnetization components (above) and the energy gap (below) versus the intensity H of a magnetic field applied in plane 
along the hard direction (X), for fixed temperature T /Tc = 0.5. Figs. 2a, b display the calculated equilibrium angle <f) (thin 
full line, referred to the right vertical axis), magnetization modulus M = (S z ) (dotted line) and magnetization components, 
(S z ) = Mcos(j> and (S x ) = Msm<j> (thick full lines), versus H. Figs. 2c, d display the calculated equilibrium angle <f> (thin 
full line) and energy gap Eq> (thick full line) versus H . The Hamiltonian parameters used for the calculations are J = 500 K, 
K2 = 5 K, 4nwf w = 1 K, S = 1, providing Tc = 670 K in mean field theory and Tc = 502.4 K in Green's function theory. 



where (S z ) = SC(x) and x = x easy = -^{g^sH 2 + 2K 2 {S Z ) + 4J(5 Z )). It turns out that uj easy monotonously 
increases as H is increased (not shown in Fig. 2c). 



For field along the hard axis, H = F X, the time varying magnetization is M = (xm x + ym y ) e 



(where x, y and z denote local axes), while the effective field is 
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(where X, Y, Z denote crystallographic axes)ii£ Expressing the effective field in local coordinates and solving the 
linearized Bloch torque equations, the resonance frequency is found to display a non monotonous behavior (see Fig. 2c), 
different for H x below and above the reorientation field = 21<2 ^ s l 
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where (S z ) = S£(x), x = x hard = -^ T (gfi B H x sincf) + 2K 2 (S Z ) cos 2 + 4J(5 Z )) and <f> is given by the equilibrium 
condition, Eq. J5J. It is worth noticing that the finite temperature mean field theory predicts a vanishing reso- 
nance frequency of the uniform mode just at H x as a consequence of the continuous in-plane reorientation of the 
magnetization. 



C. High temperature expansions 

The behavior of spin wave energies in the high temperature regime above Tq is made particularly interesting 
by the reduced dimensionality of the system. In two dimensional magnets, the spin system is highly susceptible 
to perturbations and can respond strongly even at temperatures above Tq. This means that a two dimensional 
ferromagnet in the paramagnetic regime can still produce large scale correlations when in a static applied field. 
Furthermore, these correlations are sensitive to the orientation of the applied field relative to the symmetry axis of 
the anisotropy. 

For temperatures above Tc a finite magnetization is possible only in the presence of an applied field. Furthermore, 
since the reorientation field is found to vanish at T — Tc, the field- induced magnetization will align parallel to the 
applied field. A simple estimate of the resonance frequencies above Tc can thus be made expanding the Langevin 
function for high temperatures (small x) C(x) ~ § — f^. For field along the easy axis (H = H Z Z) and along the hard 
axis (H = H X ~K.), respectively, the resulting resonance frequency is 
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with T±(T) = ^| ± j^^j^j 9l ksT ■ ^he frequencies have a linear dependence on field with different temperature 
dependent slopes. The slope is smaller for the field in the hard direction because in that case the anisotropy reduces 
the torque acting on the spins, thereby reducing the susceptibility with respect to the easy case. 
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III. GREEN'S FUNCTION THEORY 



Although mean field theory has proved useful to study reorientation transitions in multilayer films^i it is well known 
that it grossly underestimates the effect of single site anisotropy in Heisenberg models^ and that it leads to estimates 
of the critical temperature generally too high due to the neglect of spin correlations. As the reorientation field depends 
on anisotropy and on the temperature dependence of magnetization, one can expect that mean field theory is not 
capable of providing results that can be compared quantitatively with experiments in high quality epitaxial systems. 
A Green's function theory using the random phase approximationi^iS appears to be a better choice since it provides 
a more correct treatment of both issues, namely single site anisotropy and spin correlations. 

Recently, a number of Green's function theories were proposed by different authorsi2*i£*iliiiL±& to treat the problem 
of the field-induced transition in ultrathin ferromagnetic films. A thorough discussion of the limits of the approxima- 
tions used in Ref. 14 to treat the reorientation transition was provided by the authors of Ref. 0, who proposed a 
theory based on a generalization of the Caller. 19 decoupling valid for an arbitrary direction of the external field. 

In this work we present a Green's function approach which was developed independently from Ref. Il8i We start 
from the outset for rotation of the quantization axis for fields below the reorientation transition: i.e., we write the 
Green's functions equations of motion in the local reference frame. Next, by a careful treatment of the uniaxial 
anisotropy through a generalized Callen decoupling 19 of the higher order Green's functions, we obtain - in a very 
simple way even in the case of field applied along the hard axis - both the equilibrium condition and the frequency 
of the magnetic excitations at finite temperature T > 0. As a consequence of the generalized Callen's decoupling, 
the validity of our theory is limited to low values of the ratio r — K 2 j J between the uniaxial anisotropy and the 
exchange constant. For high values of r, it is advisable to use a more refined Green's function approach 12 where the 
single ion anisotropy terms are treated exactly by introducing higher order Green's functions and, taking advantage 
of relations between products of spin operators, a closed set of 2S e qua tions of motion is obtained for the anisotropy 
Green's functions. Also the quantum Monte Carlo approach in Ref. |2fJ is expected to be appropriate for large values 
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of the ratio r, while, in the opposite limit of small r, finite size effects strongly influence the 'critical region' around 
the reorientation transition). 

Coming back to our Green's function approach, we observe that complete agreement with Ref. is obtained 
for the equilibrium condition 21 while different expressions are found for the energies of the magnetic excitations^ 
Apparently, this discrepancy is due to the fact that in the Green's function theory of Ref. |18| some extra approximation 
was made, consisting in the treatment of the anisotropy in terms of an effective field, as far as the magnetization and 
the magnon energies are concerned. 

In order to study the reorientation transition and the field dependence of the magnon energies, we have taken 
particular care in correctly recovering the results of non interacting spin wave theory (see Appendix B) in the low 
temperature limit, when the local thermal averages tend to constant values, (S z ) — > S and ((S z ) 2 ) — ► S 2 . This limit 
is not recovered^ by the approach in Ref. In particular, we will show (see later Eq. I19|) that the well-known 
feature 24 of a zero frequency of the uniform (k = 0) magnon mode at the reorientation field, predicted also by free 
spin wave theory (see Appendix B), is correctly recovered by our Green's function approach in the T — ► limit. 

A. Equilibrium condition and energy of the magnetic excitations 

In order to calculate the equilibrium condition and the energies of magnetic excitations in the monolayer described 
by Eq. JQ) we use a Green's function formalism 13 supplemented by a random phase approximation (RPA) with a 
careful treatment of the uniaxial anisotropy through a generalized Callen decoupling^ Details of the calculation in 
the absence of magnetic dipole-dipole interactions are given in Appendix A. Hereafter we only quote the final results, 
where dipole-dipole interactions were included in an approximate way (see Appendix C for details). Taking a rotation 
procedure of the quantization axis z (see Fig. 1), we obtain the RPA equilibrium condition in the general case of 
in-plane field with crystallographic components H x and H z (see Eq. I|A7|I in Appendix A) 

2K 2 fs{Si) sin(t)cos(t) + gfi B {H z sincj)- H x cose/)) = (10) 

where the factor fg, 

fs = l-^p[S(S + l)-{S?Sf)], (II) 

is required in order to preserve the kinematic consistency of the magnetic excitations (see later) . 

It is worth noticing that in the T — > limit one has (S z ) — > S, so that the RPA equilibrium condition reduces to 
the equilibrium condition one would obtain in the framework of free spin wave theory, using the standard Holstein- 
Primakoff transformation for spin operators, by imposing the vanishing of the terms linear in a and a + in the 
boson Hamiltoniani^ The physical meaning of this equilibrium condition is the neglect of correlations between the 
longitudinal (z) and transverse spin components. In fact, we remind that in the framework of free spin wave theory 
the longitudinal (z) spin component has a boson representation in terms of zero- and two-boson operators while the 
transverse (x, y or, equivalently, ±) spin components are expressed in terms of one-boson operators (see Appendix B 
for details). Similarly, the RPA equilibrium condition was obtained by imposing the vanishing of the Green's function 
G Z j(w) (see Eq. <|A6(1 in Appendix A), which correlates the longitudinal and the transverse spin components. 

For field applied along the hard axis (H z = and H x ^ 0), the equilibrium angle is given by <f> = sin" 1 ( ^j^^gTj ) 
for H x < H x and <fi — J for H x > H x . The reorientation of the magnetization within the film plane occurs at a 
critical value of the applied field H x = Notice that, as in the mean field theory, the reorientation field is a 

function of temperature, but the RPA condition contains the additional temperature dependent correlation function 
(S Z S Z ) in the kinematic consistency factor fg. The latter factor is obtained by a careful decoupling of the higher 
order Green's functions coming from the uniaxial anisotropy^ and is required in order to preserve the kinematic 
consistency of the magnetic excitations: i.e., for S — 1/2, the uniaxial anisotropy does not contribute to the magnon 
energy (see later, Eqs. <|12I13I14I) '). In fact, for S — 1/2 one has fg = 0; for 5=1, which is the simplest non-trivial 
case, one has fg = \{S Z S Z ). 

The general form of the magnon energies in the momentum representation is, including also the magnetic dipole- 
dipole interaction (see Appendices A and C for details) 

- - ^ (12) 

where 

A k|| = 4J(5.f)(l - 7k|| ) + 2K 2 f s (S z ) cos 2 4> - K 2 f s (S z ) sin 2 cp + ^4irwf w (S z ) + g^ B (H x sin0 + H z cos0) (13) 
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FIG. 3: Green's function calculation of the energy gap for field applied along the hard axis (full lines) and the easy axis (dashed 
lines) at four different temperatures (T/Tb=0.5, 0.9, 1.05 and 1.2, respectively). With the Hamiltonian parameters used for 
the calculations, J = 500 K, K2 = 5 K, 4nwf m = 1 K, S = 1 (giving Tc = 502.4 K), a magnetic anisotropy of the spin wave 
energies is found to persist up to T /Tc = 1.2. 



S k|( = K 3 fs(Sf) sin 2 cj> + Unwf w {Sf) (14) 

In Eqs. H13I14|) . 7^.. = i[cos(fc x a) + cos(fc z a)] is the geometric factor for the simple quadratic (s.q.) lattice and 

w = ly9>1 ^ denotes the strength of the magnetic dipole-dipole interaction, while the factor f w — 1.078 comes from 
the calculation of dipolar sums for the s.q. lattice&i& Notice that the wavevector dependence of the dipolar sums 
was neglected for the sake of simplicity (see the Appendix C for details). Such an approximation for the dipolar 
interaction is certainly plausible in the presence of a sufficiently large gap in the spectrum. 

For S = 1 the relevant thermal averages in the RPA theory are (general expressions for any value of S can be found 
in Appendix A) 

1 + 2$(T) 

W ' = 1 + 3$(T) + 3$ 2 (T) ( 15 ^ 



((^) 2 } = S(S + l)-(l + 24>(T))(Sf) 



1 + 2$(T) + 2$ 2 (T) 



1 + 3$(T) + 3$ 2 (T) 

where the temperature dependence is self-consistently given by (N is the total number of spins in the s.q. lattice) 



(16) 



$(T) = — V 
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For zero field, in the T — > Tc limit, one has (S z ) — ► 0, so that E^.. vanishes with it and $(T) diverges. Thus from 
Eqs. ljT5)l and lfTo|) one has (5?) — > g^yy and (S Z S Z ) — > |, respectively. By equating the expression (|T7Jl (expanded 
for — > 0) with $(T) = jrgsx, one obtains the following expression for the Curie temperature of the S = 1 
ferromagnetic monolayer in zero field 



k B Tc 



ki. 



4J(1 - 7k,,) + §#2 + \4nwf u 



(A J(l - 7k|| ) + |X 2 + 4™/™) (4 J(l - 7k|| ) + §A' 2 ) 



(18) 



B. Discussion 



In Fig. 2b we show the self-consistently calculated equilibrium angle cj>, magnetization M — (S z ), and magnetization 
components (S z ) = (S z ) cos</> and (S ) = (S z )sm(j), as a function of the magnetic field H, applied along the hard 
axis, for fixed finite temperature (T/Tc = 0.5). Correspondingly, in Fig. 2d the self-consistently calculated energy 
gap Eq is shown as a function of H. 

In contrast with the mean field result shown in Fig. 2a, in the framework of Green's function theory we find that the 
orientation angle <p varies discontinuously at the reorientation field and that the magnetization modulus M = (S z ) 
decreases as H is increased from to the critical value H x = 2K2 ^ s ^ Si I . This is related to the more accurate treatment 
of spin correlations in the Green's function theory with respect to mean field theory. In fact we observe that both the 
temperature and field dependence of (S z ) occurs through the quantity $(T) in Eq. (|17|) . which is determined by a 
summation over all magnon modes. At fixed temperature, the strongest contribution to $(T) comes from the k = 
mode since this has the lowest energy for H x — > (H x )~. This explains why M = (S z ) decreases in this limit. In 
turn, the orientation angle 4>, whose value is determined by the value of (S z ), undergoes an abrupt variation at the 
reorientation field. 

As a consequence of the discontinuity in the orientation angle <f>, the energy gap Eq displays a discontinuous jump 
at the reorientation field H x = 2K2 g^ Si ' ( see Fig. 2d) and the energy does not go to zero at the reorientation field. 
This is in contrast with the finite temperature mean field result displayed in Fig. 2b. 

It is important to note that the softening of the k = magnon energy at the reorientation field 2 ^ is correctly recovered 
by our Green's function approach in the T = limit. In fact for T — > and S = 1 one has (S z ) — (S Z S Z ) — > 1 and 
the energy of the k = magnon mode exactly coincides with the result of non interacting spin wave theory for field 
H x applied along the hard direction (see Appendix B for details) 



E = [{K 2 f - (gfi B H x ) 2 ] 



for H x < H? 



K 2 

E = [g^ B H x - K 2 ] h x [g^ B H x + Anwf w ] * , for H x > H x . (19) 

This energy is immediately seen to vanish continuously at the reorientation field H x = . 

The magnitude of the gap energy at zero field is interesting to consider. The zero field gap value is a measure of 
the effective anisotropy energy, and is strongly temperature dependent: the zero field gap value vanishes for T > Tc, 
as well as the reorientation field. This is clearly a consequence of the Callen decoupling^ of the higher order Green's 
functions generated by the uniaxial anisotropy, yielding a decoupling factor proportional to the magnetization, see 
Eq. I|A14(1 . Physically this is related to the absence of a collectively ordered magnetic state for T > Tc- However, 
as discussed by Jensen et al.£ the vanishing of the effective anisotropy for T > Tc does not indicate that the 
microscopic anisotropy vanishes either: a single magnetic moment in the paramagnetic state is still subject to the 
anisotropy even if the net magnetization is zero. While in Ref. such a "paramagnetic anisotropy" was revealed 
by magnetic susceptibility measurements, in the present paper we rather focus on the gap energy behavior. This 
feature is illustrated in Fig. 3, where the gap energy is plotted as a function of the magnetic field applied in the 
easy and in the hard direction, at various fixed temperatures, both below and above Tc- For T > Tc, the easy- and 
hard-direction energies are each zero at zero field, and their field dependence becomes more and more linear as T 
increases. However, the calculated slopes are different, in agreement with the mean field argument: the hard axis 
alignment slope is smaller than the easy axis alignment slope. Finally, it is worth noticing that for a realistic value of 
the ratio K 2 /J in an ultrathin film (w as in our numerical calculations), the magnetic anisotropy of the energy 
gap is expected to persist well above Tc (up to T/Tc ~ 1.2, see Fig. 3d). 
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IV. CONCLUSIONS 

Using a Green's function theory with random field approximation we have shown that the field driven reorientation 
transition of a ferromagnetic monolayer with uniaxial anisotropy is expected to be continuous at T = (in agreement 
with free spin wave theory) and discontinuous at finite temperature. In contrast, mean field theory predicts a field 
driven reorientation transition continuous at any temperature. Below the critical temperature Tc, the Green's function 
theory predicts a magnetic excitation gap (the resonance mode) which does not completely soften at the reorientation 
field. 

The question of the order of the transition should be considered as a still open question, because it is not fully clear 
if and how the decoupling scheme affects that behavior. One step in this direction might be accomplished performing 
a decoupling of the Green's functions on a higher level. 25 It would also be interesting to analyze the order of this 
type of transition in experimental systems (e.g., by ferromagnetic resonance measurements of the uniform mode as a 
function of the field). 

The magnetization of a two dimensional ferromagnet is highly susceptible above the critical temperature Tc. This 
means that a large magnetization can be induced with relatively small applied fields in the paramagnetic regime. Such 
effects were demonstrated in beautiful experiments on ultrathin Fe/W(110) films by Back et ali More recently, such 
"paramagnetic anisotropy" effects were investigated by means of magnetic susceptibility measurements on ultrathin 
Co/Cu(l 1 17) films by Jensen et al. 2 In our case, this means that a two dimensional film with well defined anisotropies 
can support magnetic excitations in the paramagnetic regime with energies that depend upon the orientation of the 
applied field relative to the anisotropy symmetry axis. We have demonstrated this with our Green's function theory 
in the RPA and shown that there can be a sizeable difference in magnetic excitation energies at finite fields. We show 
that the field dependence of the energy is strongly temperature dependent and different for field alignments along 
easy and hard directions. 

It would be interesting to put our prediction of an anisotropy of the magnetic excitation energies in the paramagnetic 
phase to an experimental test. This may be possible in ultrathin magnetic films using inelastic light scattering or 
ferromagnetic resonance measurements. 

In fact, recent Brillouin light scattering experiments in Fe/GaAs(100) ultrathin films (tp e = 4.5 A) showed 
anisotropy in observed frequencies above some ordering temperature^ However, it is not clear from the results 
whether the film is continuous or is instead a collection of superparamagnetic islands. 

Our mean field arguments, for example, apply equally well to paramagnetic films above the Curie temperature and 
to collections of superparamagnetic particles above the blocking temperature. The only difference might be that a 
superparamagnetic "film" would be comprised of a distribution of particle sizes, each with a different local internal 
magnetic field. Resonance in the superparamagnetic system would contain many different frequency contributions 
and should therefore have a correspondingly broad linewidth. Resonance in the paramagnetic continuous film would 
have a linewidth determined by higher moments in the spin correlation functions, and should therefore also have a 
distinctive temperature dependence different from the superparamagnetic case. 
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APPENDIX A: GREEN'S FUNCTION FORMALISM 

In order to calculate the equilibrium configuration of the monolayer magnetization and the energies of the magnetic 
excitations, we consider the following Green's functions 

Gg. = ((Sf;Sr-)>; a = +,-,z (Al) 

where S t = Sf ± iSf and x,y,z denote a reference frame where z is the equilibrium direction of the in-plane 
magnetization, forming an angle cf> with the easy in-plane direction Z (see Fig. 1). The equations of motion for the 
time Fourier-transformed Green's functions Gf (lu) are 



u,G«H = ([Sf, SJ]) + «[$?,«]; Sr» 



(A2) 
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where (A) = i?Tr(e 131-1 A) denotes a thermal average (/3 — -jrKp and Z = Trie l37i ) is the partition function). Taking 
into account that S x — S x cos <p + S z sin <j> and S z = —S x sin <p + S z cos <fr, one obtains 

lo Gf.(w) = ([Sf,Sr]) 

=F J y~]{{S*Sk — Sp,Sf\S- )) u 

k 

± K 2 cos 2 ^{{SfSf + S±S*; Sr)) u 

T ^2 sin 2 5r + firfi* + S?S+ + S7)) w 
± K 2 sin </> cos <j> 

((2S*Si — S^S^ — -{S^S i + S i Si~); Sj )} u 
± ffMs^ Z (sin G*j (w) + cos G± (w)) 

± ff/ is^ X (sin0G±H-cos0G?.H) (A3) 



w G?.(o;) = <[S?,S7]) 

1 k 

\k 2 sin 2 0((S+S+ - S-S-- SJ)) U 

+ 7^2 sin cos 

{{SfSf + St $t — Si S* — S z Si ;Sj 
+ \ gilB H z sin0(G*(w) - Gy( w )) 

- ^Bff^cos^GtH-Gr.^)) (A4) 

The higher order Green's functions on the right hand sides of the previous equations have now to be decoupled 
to obtain a closed set of equations. We start using the naive RPA decoupling (a, = +,—,z) for all the Green's 
functions (i.e., both for i ^ k and i — k; kinematic consistency corrections will be taken into account later on) 

((S? S f l ■ SJ)) U « (Sf }Gf » + (5f)G« ( W ) (A5) 

Next, taking into account that the local axis z is, by its own definition, the equilibrium direction of the magnetization, 
we can safely put (S + ) w and (S~) w in the equations of motion of the three Green's functions G^(w), G~-(w) 
and G?(a>). 

The equation of motion for Gf -(w) provides the equilibrium configuration. In fact in the previous approximations 
we obtain 

<* G ii ( w ) ~ ^ [ G 5 M - (w)] x [ tf 2 (Sf } sin(20) + g m (H z sin - H x cos 0)] ( A6) 

so that the equilibrium angle cj> is obtained imposing the vanishing of the term in braces on the r.h.s. 

2K 2 {S*) sin <j> cos <f) + gii B {H z sin <f) - H x cos <j)) = (A7) 

For field applied along the easy axis (H z ^ and H x = 0), one obviously obtains = 0: i.e., the magnetization 
is always directed along the field (Z = z) direction. For field applied along the hard axis (H z = and H x ^ 0), the 
equilibrium angle is given by 

* = 8i "~'<itS>'- fo '- ffj ^ 

4, = |, for H x > H; Y (A8) 

so that a reorientation of the magnetization within the film plane occurs at a critical value of the applied field 
H x = 2K g^ - which turns out to be temperature dependent. For H x > H x , the magnetization becomes parallel to 
the external field H x applied along the hard direction. 
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The equations of motion for GfAio) and G7(w) provide the frequencies of the excitations with respect to the equi- 
librium configuration. In fact for the space Fourier-transformed Green's functions G ± (lo, ky) = G^ (u;)e~ 4k ii ' riJ 
(where ky = (k x ,k z )) we obtain 

±u(s?)(i- lkll )G ± (u,yi ll ) 

±2K 2 (S*)cos 2 ^(w.kn) 

TK 2 (Sf) sin 2 [G+(lo, k,| ) + G-{u, ky)] 

±^i B [i7 z cos</>-f-i7*sin</>] G ± (w,ky) (A9) 

where J 7 ' ± (ky) = e _lk "' riJ ([5 f ,57]) and ^ = ^— J2s e~ lkl, ' S ($ are the vectors joining a given lattice site to 

its z nn nearest neighbors). For a square lattice with lattice constant a one has z nn = 4 and the geometrical factor 
takes the form 7^.. = ^[cos(k x a) -I- cos(fc z a)]. 

The previous equations can be rewritten in a more compact form 







"G+(w,ky) " 




> + (k||) " 


-Bk„ w + A k|| 




G-(w,ky) 




.^"(ky) . 



(A10) 



where the quantities Ak y and Sky take the following expressions 

A k|| = 4J(5?)(1 - 7k„ ) + 2K 2 (Sf) cos 2 - # 2 (5?) sin 2 tf> + gfi B (H z cos</» + H x sin0) 

Sk„ = ^2(5?) sin 2 (All) 

The energies of the magnetic excitations with respect to the equilibrium state are then given by 

El, =(^ k|l ) 2 =A k|| -B^ (A12) 

It is worthwhile noticing that in the zero temperature limit one has (Sf) — > 5 so that the free spin wave frequency 
is correctly recovered. In particular, in the T = limit we are able to recover the well-known feature of a zero 
frequency mode^ for field applied along the hard axis and equal to the critical reorientation field, in contrast with 
the results obtained by the Green's function approach in Ref. It is worth noticing that other Green's function 
approachesi^i^i^ were not able to recover such a feature, either. 

So far, in our calculations, the problem of the kinematic consistency of the magnetic excitations was neglected 
since the naive RPA decoupling was performed also for the Green's functions coming from the anisotropy term in 
the Hamiltonian. In the case of field applied along the easy axis (H z ^ 0, H x = 0), it is well known that a more 
correct treatment of the anisotropy term is obtained by using the Anderson-Callen decoupling to decouple the equal 
site [i = k) Green's functions^ 

«5±5? + 5?5f ; 5r)) w « 2<Sf)G±(u;) x jl - ^ i S ( S + (W)]} (A13) 

In this way, for 5 = 1/2, the uniaxial anisotropy does not contribute to the frequency of the excitations, as expected, 
since for 5=1/2 the anisotropy term in the Hamiltonian is a constant. In the case of field applied along the hard 
axis (H z = 0, H x 7^ 0) or in the general case (H z ^ and H x ^ 0), we assume the same prescription for kinematic 
consistency to be valid as in the easy case, i.e. we replace 

K 2 (Sf) - K 2 (Sf) [l-g^a W S + !) " ( S i S i)}} (A14) 

both in the equation (|A7() providing the equilibrium condition and in the equations IjAllI) and (1 A 12(1 giving the 
energies of the magnetic excitations. In the case 5=1 considered in this paper, the prescription for kinematic 
consistency reads simply 

2K 2 (Sf)^K 2 (Sf)(SfSf) (A15) 

In order to determine the average (5f Sf), we follow the method exposed by Calient in his paper. In addition to the 
Green's functions Gjy = ((Sf; 57)), we consider the Green's functions Qf^ a = ((e aS * Sf; 57)) where a is a parameter. 
The (?'s are found to satisfy the equations of motion 



W-A(ky) B(ky) 




'0+(w,k,|) " 




'^(k||) ' 


-5(ky) W + A(ky) 




.Sa(^ky) _ 




.^„-(ky) . 



(A16) 
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which differ from the equation of motion for the G"s only for the terms on the r.h.s., defined as J-^(k^ 
-L\\. 



o— tkii 



{[Sf, e a * Sj ]). Thus from the spectral theorem 



IB A) = — lim 

v ' 2tt 



+00 j 

-3— S) ) ( W + ie, k,| ) - ((A; B) ) {u> - it, k„ )] 



(A17) 



applied to G + and one can obtain, respectively, the thermal averages 

(Sr S +) = F+(k H ) $(T) 
(e aS *5r5+> = ^+(k|) $(T) 



where 



$(T) = — V 



k N L 



^k„ 1 
-^coth(-/?£ k|| )-l 



For the two-spin correlation function we obtain 



(5" 5+) = 2(S*MT) 



(A18) 



(A19) 



(A20) 



Following CallenjiS it is now possible to find a second-order differential equation for .F+(kii) as a function of the 
parameter a that, when solved with opportune boundary conditions, allows to obtain the magnetization as a function 
of $(T) 



[5 - *(T)][1 + $(T)] 2S+1 + [5 + 1 + $(T)][$(T)] 2S+1 



and finally 



[1 + $(T)] 2S +! - [$(T)] 2S +! 
(5*5?) = 5(5 + 1) - (5r5+) - (5?) = 5(5 + !)-(! + 2$(T))(5?) 



(A21) 



(A22) 



APPENDIX B: FREE SPIN WAVE THEORY 

We think it useful to briefly present here the results of non interacting spin wave theory, both for the equilibrium 
condition and the energy of the magnetic excitations in an anisotropic monolayer described by the microscopic 
Hamiltonian Q, since in the T — » limit the spin wave results are expected to be exactly reproduced by Green's 
function theory. Denoting by x,y,z the local reference frame where z is the equilibrium direction of the in-plane 
magnetization, forming an angle 4> with the easy in-plane direction Z (see Fig. 1), the local spin components are 
expressed in terms of boson operators via the Holstein-Primakoff transformation^ 



Si 




afaj 
25 



S- 




5f 



5 — OjCLj 



(Bl) 



while the spin components in the crystallographic reference frame are obtained from the local ones through the 
equations S x = S x cos<j> + 5 z sin</> and S z — —S x s'm(f> + 5 z cos0. In the following, for the sake of simplicity, we 
neglect both dipolar interactions and kinematic consistency corrections to the anisotropy term in Eq. Q . Substituting 
in £Q| , expanding the square roots in Eq. IjBljl (such an approximation is valid for low values of the occupation number 
rij = a^aj, that is for low temperatures), and keeping only up to quadratic terms in the boson operators one obtains, 
after a space Fourier transformation aj = — 7= a^,,e~ 



•ikll -Ti 



H ~ Hi +H 2 = \/|-X^ ak ii +a k,|) [2if 2 5sin ( />cos(/> + .g^ B (ir 7 sm ( />- H x cos0)] 



E 



-Aic.o,^ awn H — Bit,, (die,, ct w ,, -h tit ct~*~i ) 

Km ^ K 1 O K V K K Km — JC II / 



(B2) 
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where 

A k|| = 4J5(l-7k„) + 2K 2 S cos 2 (/}- K 2 S sin 2 <p + g^i B (H x sine/) + H z cos <p) 

S k ,| = -K 2 S sin 2 0. (B3) 

The equilibrium condition is obtained by imposing the vanishing of the one-boson Hamiltonian Tii 

2K 2 S sin<pcos(/) + gn B (H z sine/)- H x cose/)) = (B4) 

It is readily seen that the spin wave theory result Eq. I|B4(I is reproduced by Green's function theory in the limit 
T -> 0, see Eq. (|X7|) . 

The two-boson Hamiltonian Ti 2 can be readily diagonalized£i& and the energy of the non interacting spin wave 
excitations turns out to be 

^i, = - (B6) 

Taking into account both dipolar interactionaSi^ and kinematic consistency corrections^ for field H x applied along 
the hard in-plane direction one finally obtains the free spin wave energy explicitly reported in Eq. (|19fl . 



APPENDIX C: CONTRIBUTION OF DIPOLAR INTERACTIONS 



The inclusion of magnetic dipole-dipole interactions (i.e., the last term in Eq. in the equations of motion for the 
Green's functions leads to a considerable amount of additional calculations. As they are tedious but straightforward, 
hereafter we only quote the final result after RPA decoupling. First we consider the commutator of Tidip with Sf in 

2 2 

order to investigate the contribution of the dipolar interaction to the equilibrium condition. We have (w = 9 3 ) 

3 



XX 
kl 



d zz ) sin (j> cos < 



(cos 2 <j) 



G+.H-G--H 



XY 
kl 



sm ( 



,YZ 

1 kl 



COS < 



G+M + Gr.( w ). 



2i 



where the dipolar sums are defined as (a, /3 = X, Y, Z) 



D 



a0 



(kn)-E d : 



a/3 ik|| 
ij e 



IT 



a 
r ij r ij »k|, 
2 

r ■ 



(CI) 



(C2) 



For a monolayer with the crystallographic Y axis normal to the film plane one has rL = so that D = D = 0; 
these relations are approximately satisfied even in an ultrathin film. Moreover, in a lattice with a center of inversion 
symmetry one has D xx = D zz and D xz = 0. Thus, as expected, it is proved that the dipolar interaction does not 
contribute to the equilibrium condition of an in-plane magnetized ultrathin film. 

Next, we consider the commutator of Tidip with Sf 1 in order to obtain the contribution of the dipolar interaction 
to the magnon energies. In the same approximations as before we obtain, after a Fourier transformation 



j^dip 



1 



X 



1 - 3 



sm i 



rfj cos < 



1 - 3 



' r x cos < 



r z sin < 



> r x cos ( 



r z sin < 



(C3) 



In order to simplify the numerical calculations, we neglect the dependence of A^ p and on k|| . Taking into account 
that for a two-dimensional lattice with a center of inversion symmetry one has 



e4 ( rl ) 2 =y:i C 1 ) 2 =- 1 -y:i=--«f« 



(C4) 
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with f w = 1.078 for the s.q. lattice £tL& we thus obtain the approximate expressions i|13|) and l|14fl reported in Section 
II. Clearly, such an approximation on the dipolar part of the Hamiltonian is allowed as long as the strength of the 
dipolar interaction w is smaller than the intensity of the uniaxial anisotropy and/or the Zeeman term. 
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